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Abstract 



Q I We simulate the low temperature behaviour of an elastic chain in a random po- 



tential where the displacements u{x) are confined to the longitudinal direction (n(x) 



^ ' parallel to x) as in a one dimensional charge density wave-type problem. We calculate 

the displacement correlation function g{x) = {{u{x) — u(0))^) and the size dependent 
average square displacement W{L) = {{u{x) — n)^). We find that g{x) ~ x^^ with 
r] ~ 3/4 at short distances and t] ~ 3/5 at intermediate distances. We cannot resolve 
the asymptotic long distance dependence of g upon x. For the system sizes considered 
we find g{L/2) ocW^ L^^ with x — 2/3. The exponent r] ~ 3/5 is in agreement with 
the Random Manifold exponent obtained from replica calculations and the exponent 
X — 2/3 is consistent with an exact solution for the chain with transverse displacements 
(u(x) perpendicular to x). The distribution of nearest distances between pinning wells 
and chain-particles is found to develop forbidden regions. 

PACS numbers: 74.60.Ge, 61.41. +e, 68.10,-m 



1 Introduction 

An elastic medium in a random static potential models numerous physical phenomena such 
as domain walls in impure magnets, ^, flux line lattices in dirty Type II superconductors 
in the mixed state, 0, ^ polymers in a disordered environment, etc. The configuration 
of the elastic medium is determined through the competition between the elastic energy of 
the medium and the pinning energy due to the inhomogeneities. This is a confiict between 
ordering and disordering. In less than four dimensions the random potential will win leading 
to loss of translational order at zero temperature. 

The degree of distortion of the elastic structure grows with the considered length 
scale. Two characteristics are often used to quantify the distortion. One is the displacement 
correlation function 

^(r) = (Kr + ro)-n(ro)f) (1) 

where u{r) denotes the displacement -u of a point in the elastic medium from its ideal lattice 
position r and the angular brackets denote a thermal average, average over choice of origin. 
To, and an average over the disorder. 

This function is used in collective pinning theory to calculate the correlated volume, 
Vc [^, as well as in the theory of collective creep [0]. Furthermore, it has been used in 
attempts to calculate the correlation function of the translational order of the elastic system 



There have been many attempts to calculate the displacement function either by 
use of scaling arguments (see for instance Ref. |^, |10|, 0) or using replica techniques |]11, 



T2| , |T3], [T^. These studies obtain three separate regimes for g{r). The first regime is the 
Larkin regime and is defined by g{r) < R^, where Rp is the length scale of the pinning 
potential. This short distance regime describes the behaviour within a correlated volume 
and is expected to be derivable from perturbation theory 0, [^ . The second regime is often 



referred to as the Random Manifold regime, which is defined over the range R^ < g{r) < Oq, 
where oq is the ideal lattice spacing. The final regime is the asymptotic long distance 
regime, defined by g{r) > Oq. The displacement function scales differently in each of the 
three regimes. 

There are three important dimensionalities in the problem. Firstly, there is the 
actual dimension of the space that the random potential occupies, D. Secondly there is the 
dimension of the components of the elastic medium, d (i.e. points are d = 0, lines d = 1, 
planes d = 2), and finally there is the number of dimensions perpendicular to the components 
of the elastic medium that the displacements move in, n. For example, D = 2, d = n=l 
represents an interface in a two dimensional random bond Ising model, whilst D = 3, d = 2, 
n = 1 is a planar interface in three dimensions, and D = 3, d=l, n = 2 represents flux 
lines in dirty type II superconductors. In general, D = d + n. The precise scaling in each of 
the three regimes of g{r) will usually depend on all three of these dimensionalities. 

Another measure of the length scale dependence of the distortion is obtained from 
the fluctuations in the displacement field averaged over the whole system and then studied 
as function of different system sizes L. I.e., 

W{L) = {{u{r) - {u)r) (2) 

An exact solution for W{L) was offered by Huse, Henley, and Fisher |l|. This was 
for an interface in the continuum version of the two dimensional random bond Ising model 
(i.e. D = 2,d = l,n = 1). They derived an exact scaling exponent for the size dependence 
of W{L). Namely W{L) ~ L^^ with x = 2/3. The discrete version of this model is similar 
to ours (discussed below) in that both consist of a chain of particles but the above model 
is subject to transverse displacements whilst the displacements in our chain are confined to 
the longitudinal direction. 

Here we consider a one dimensional chain of particles connected by elastic springs 



embedded in a static random potential (i.e. D = l,d = 0,n = 1). Periodic as well as open 
boundary conditions are studied. We calculate the displacement function in Eq. |I| and the 
average square displacement in Eq. ^by Monte Carlo annealing. We find that g{x) exhibits 
huge fiuctuations and that the averaged quantity is very sensitive to the specific details of 
the annealing and averaging procedure used. Hence, one has to be very careful in order to 
make sure that thermal equilibrium is achieved and that the results are reproducible. 

Our motivation is to compare our results to the replica calculations of Bouchaud, 
Mezard and Yedidia (BMY).[|Tl|, [T2| For the given dimensionalities, they obtain the following. 



In the Larkin regime, they obtain g{x) ~ x^^ with rj = 3/2. |T^ In the Random Manifold 
regime, they get the same scaling form as the Larkin regime but with an exponent rj = 3/5. 
Finally, in the long distance regime they obtain g{x) ~ x^^/i/loga; with r] = 3/4. 

In the case of periodic boundary conditions we observe two regimes of different growth 
of g{x). We can fit our data to an algebraic form g{x) ~ x'^^'. In the short distance regime 
rji ~ 3/4, in disagreement with the calculation of BMY, and in the intermediate distance 
regime 772 — 3/5 in agreement with BMY's replica calculation. ||TI|, |12| We cannot comment 



on the scaling in the long distance regime since computational restrictions prevent us from 
considering long enough chains. 

In contrast to the theoretical expectation we find that the exponents measured in 
the simulation do depend on the strength of the random potential, with 772 approaching 1/2 
for the stronger potential. This could possibly be a finite size effect. 

It is very difficult to produce satisfactory statistics for the open chain. However, our 
simulation indicates that there is no difference between the open and the periodic chain. 
One might anticipate a difference, since in one case the density of particles in the chain 
is allowed to change whereas in the other case the density is fixed. This difference will of 
course be most pronounced at long distances when the relative displacement becomes of 
order qq. Hence one could argue that our chains are too short to resolve the difference. Our 



largest relative displacement is of order 2ao. 

We also find that our simulations give W{L) ~ L^^, with x — 2/3. This is consistent 
with the exact solution for the transverse chain. One might expect the two systems to 
have different exponents, but it is worth noting that for disorder free cases where the dis- 
placements are due to thermal effects only, the two systems behave with the same thermal 
roughening exponent oi t] = 1/2. It is also worth noting that super-universal exponents can 
occur. Numerical simulations of Kardar and Zhang show that the exponent x = 2/3 holds 
for not only the exact solution where D = 2,d = l,n = 1, but also for D = 3,d = l,n = 2 
and D = 4,d = l,n = 3. Finally, the difference between transverse and longitudinal dis- 
placements will be most important for relative displacements large compared to the lattice 
spacing. We can not exclude a cross-over to a different x-exponent at much longer distances. 

In order to quantify the amount of correlation between the particles in the elastic 
chain and the randomly positioned pinning wells we have studied the distribution of nearest 
distances, P(r), between particles and pinning centres. We observe that P{r) develops gaps 
around values of the order of the range of the pinning centres. These forbidden regions vanish 
with increasing temperature. This behaviour can be understood in terms of an effective one 
particle model. 

In addition to the spatial behaviour of the system we also investigate the relaxation 
times of the system. We have looked for the possible existence of ergodicity breaking by 
measuring the relaxation times as function of system size as usually done in spin glasses 



1% . No ergodicity breaking was observed. 



2 The Model 

Our model consists of a one dimensional chain of A^ harmonically coupled particles with 
positions Xi. The equilibrium ideal lattice spacing is qq and the coupling constant of the 



elastic interaction is k. The chain is then subjected to a static background potential created 
by Np attractive pinning wells with random positions xf*". The resulting Haniiltonian for 
the system is 

N 1 N 

H = J2 7^{x^+l - a;. - ao)' + E V{x^) (3) 

We use oo as our length unit and k as our energy unit, thus setting ao = k = 1. The random 
potential at a point x is V{x) and it is generated by 

n^) = -E/.exp M^ ^^ M (4) 

The pinning wells are represented by attractive Gaussians of strengths Ap and ranges Rp. 
We always take the density of pinning centres as Up = Np/N = 1. 

We perform Monte Carlo annealing simulations on this Hamiltonian in order to study 
the low temperature properties of the system. 

We now describe the simulation. Standard Metropolis importance sampling is used 
p9|] . The initial configuration for the chain is chosen at random. We then choose a particle 
at random and give it a random displacement of —7 < dx < '-/. The factor 7 depends on 
the temperature and is chosen such that the acceptance rate is kept at about 60% for all 
temperatures. When calculating g{x) we start from a high initial temperature Tj and anneal 
the system down to the final temperature Tj. We use an inversely linear time dependence 



of the temperature [^ 



r(t) = ^ for t>0. (5) 

The time unit in the above equation is given by ngweep Monte Carlo sweeps through the 
lattice. Thus the time unit At = 1 corresponds to Nrigy^eep Monte Carlo steps. We find that 
nsweep = 400 is sufficient since the results are the same for larger values of Usweep- 

At the end of the cooling we check that the simulation a.t T = Tf is in equilibrium by 
calculating the specific heat C = {{H'^) — {H)'^)/kBT'^ (here ks is Boltzmann's constant). 
At the very low final temperature we expect C to be given by the equipartition value 
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C = NkB/2. We then perform a thermal average of g{x) by samphng the configuration of 
the particles every few sweeps through the lattice. 

The results for g{x) is supplemented by the distribution, P{r), of distances r between 
the particles in the chain and the pinning centres nearest by. This distribution contains 
directly information about the degree of correlation between the particle position and the 
positions of the pinning centres. 

In addition to the calculation oi g{x) we also measure relaxation times in simulations 
performed at fixed temperature. The relaxation times are obtained from the evolution of 
the elastic part of the energy, E^i/N, in Eq. ^. In these simulations we start out from a 
random lattice configuration. We then measure the elastic energy per particle as function 
of Monte Carlo time, E^i{t)/N [0. 



3 Results and Discussion 

Firstly, we study a system of length 200 with Ap = 0.01 and Rp = 0.1. The displacement 
function, g{x), scales differently depending on whether the fluctuations in the displacement 
field are due to thermal roughening or disorder roughening. For a finite temperature, Tf, 
there is a temperature dependent crossover length. It Ijl^, that separates the two scaling 
regimes. For x < It thermal effects dominate and g{x) goes linearly in x. For x > It 
disorder induced displacements dominate. Thus for any finite temperature, disorder will 
always dominate if long enough length scales are considered. For the considered chain, if 
we consider final temperatures Tf ^ Ap, the pinning potential becomes irrelevant since the 
length scale It is greater than the system size and g{x) is given by the exact expression 
for the case where there is no disorder in the system and the only fluctuations are due to 
temperature: 

g{x) = kBTfx{l - ^) (6) 



where N = 200 is the length of the chain. For Tf <^ Ap, the pinning potential totally 
dominates, since It is less than the lattice spacing, and g{x) deviates from the expression 
in Eq. ^. If we then increase the temperature, we can clearly identify a change in the short 
distance behaviour towards the linear behaviour of the pure system. Unfortunately large 
fluctuations limit the systems sizes affordable. For this reason we have not been able to 
resolve the temperature dependence of It- Below we concentrate on the low temperature 
region T -C Ap. 

We show in Fig. 1 g{x) measured on a periodic system of length A^ = 200 with 
Ap = 0.01 and Rp = 0.1. The solid curve is obtained by averaging over 600 different 
realisations of the random potential using the same annealing procedure for each realisation. 
(The dotted line is discussed below). In order for the Monte Carlo method to be able to 
explore as much of the phase space of the Hamiltonian as possible, we choose an initial 
temperature Tj ^ Ap so that the disorder is irrelevant and thermal effects totally dominate 
in the initial stages of the simulation. The initial temperature here is Tj = 1. The solid 
line corresponds to a final temperature of Tf = 10~^. We checked that an average over 800 
realisations gives the same curve where as 400 relaisations were significantly different. A 
huge number of realisations are needed in order to obtain satisfactory statistics. This is 
due to the low dimensionality of the problem. Fluctuations are very significant as is seen 
from Fig. 2 where we have plotted a host of the g{x) curves for a randomly selected subset 
of realisations of the random potential. These strong fluctuations are already present in 
the pure system (no random background potential). We studied the pure system in order 
to check that our annealing procedure correctly reproduces the exact result in Eq. ^ for 
Ap = 0. 

We also addressed the question of self-averaging. We did this by comparing the result 
obtained for g{x) by averaging over many different realisations of the disorder to the result 
obtained for one single realisation of the disorder averaged over the initial configuration. 
The result for the average over initial configurations is shown as the dotted line in Fig. 
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1. The final temperature is Tf = 10^^. The two different types of average give the same 
qualitative results. However, we attribute the slight difference in the curves to a finite 
size effect. Since the curve averaged over initial configurations is for a single realisation 
of 200 pinning centres, we can expect that averages of this nature will not sample the full 
distribution of pinning configurations as effectively as the same average for a larger system 
or for averages over different realisations of the disorder. Thus, we expect that the two types 
of averages will converge in the thermodynamic limit. I.e., the chain is self-averaging. 

In Fig. 3 we show a double logarithmic plot of g{x) for a periodic chain of length 
A^ = 10^ with Tf = 10"'^ for different values of Ap and Rp. We also show the result for an 
open chain. One can consistently identify two algebraic regimes where we define x = Xc as 
the crossover distance and g{x) ~ x"^^^ in both regimes. We compare these two regimes to 
the Larkin regime and the Random Manifold regime of Ref. |ll||. In Ref. [^ they obtain a 



third asymptotic regime for which g{x) > a^, where they obtain g{x) ~ x^^^/y/\ogx. It can 
be seen from Fig. 3 that the region g{x) > a^ is too small to obtain any information about 
whether the predicted scaling in this regime is correct or not. To study this we would need 
to consider much longer chains which we are prevented from doing due to computational 
limitations. The short distance Larkin regime a; < Xc is in all four cases well described by an 
exponent rji ^ 0.73. This differs from the value exponent ?7 = | expected by BMY. [|^. The 
exponent r/2 in the Random Manifold regime x > Xc is about 0.63 for Ap = 0.01. See the 
solid line (-Rp = 0.1) and the dashed line (-Rp = 0.2). The arrow marks the value of Xc for 
Ap = 0.01 and Rp = 0.1. It should be noted that the crossover distance Xc always is found 
to be only a few lattice spacings. The deviation of this exponent from analytic expectations 
might be due to the inapplicability of continiuum elasticity theory at such short length 
scales. When Ap is increased to 0.05 the exponent 772 decreases to about 0.57 as is seen from 
the dot-dashed curve. These two exponents are very close the exponent r^ = | as predicted 



in Ref. [11 1. The open chain appears to behave identically to the periodic chain. See the 



dotted curve which describes an open chain with Ap = 0.01 and -Rp = 0.1. The crossover 



distance Xc is predicted to be dependent on the pinning parameters Ap and Rp. The figure 
illustrates the difficulties one has in identifying any functional dependence of Xc{Ap,Rp). 
It is also assumed [^ |^, |ll| that g{x) at the crossover distance should scale as g{xc) ~ R^- 
Again, it is difficult to extract concrete information about this. 

It is interesting to point out that an exponent 772 = 1/2 is obtained perturbatively by 



considering a continuum Hamiltonian with a delta-correlated random potential. |T^ Consider 

H = j dx{]^K[^^f + V{x + u{x))}. (7) 

The extrema condition SH/Su = takes the simple form 

= -.^ + V'ix + u{x)) ^ -.^ + V'ix) (8) 

The approximation is correct to lowest order in V. By integration we obtain 

u{x) = - r dx'V{x') + u{0). (9) 



K Jo 



The displacement function in Eq. |l| is then given by the following pretty expression 

g{x) = — r dxi r dx2{V{xi)V{x2)). (10) 



K^ Jo Jo 



Hence, if the random potential is delta correlated g{x) ~ x^^ with 77 = 1/2. We measured 
the correlation function of the background potential in the configurations generated by the 
simulation and found that the correlations decayed to zero over a very small number of 
springs in the chain. This decay distance was of about the same value as the crossover 
distance of g{x) (see Fig. 1). We were unfortunately not able to identify the detailed 
functional form of (y {xi)V {X2)) ■ This argument suggests that the exponent r/i k, 0.73 
measured in the short distance regime is an effect of the short distance behaviour of the 
effective correlations of the potential. This perturbative result might also have some bearing 
on the behaviour observed in the intermediate regime. At these distances the random 
potential can be considered delta-correlated and if Eq. |10| were applicable one should expect 
a value for the 77 exponent of 1/2. This perturbative value is not very different from the 
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actually measured 772 ~ 0.57 to 0.63. Nor is it very far from the replica prediction rj = 3/5. 
Perhaps it is appropriate to mention that perturbation theory is usually only assumed to 
be applicable as long as g{x) < R^. This is probably correct when the random potential 
consists of a very high density {Np ^ A^) of very sharp pinning centres. In the case we have 
simulated Np = N. Only few particles are positioned in the steep regions of the pinning 
wells where a displacement of order Rp produces significant changes in V{x + u{x)). Hence, 
the approximation used in Eq. || might still be reasonable for most of the particles even 
when g{x) becomes equal to or somewhat larger than R^. 

We show in Fig. 4 a double logarithmic plot of g{L/2), together with W{L) (see Eq. 
^ for parameters Ap = 0.05, Rp = 0.1 and Up = 1. The two solid lines have both slope 
4/3. One sees that the behaviour of both quantities are consistent with the exact result 
for the transverse chain W{L) ~ L^/^.[0 However, one notices also that the data points at 
L = 10'^ fall below the straight line. This might indicate a crossover to a different behaviour 
at longer distances. Hence, we cannot exclude that the asymptotic behaviour of W and 
g{L/2) for the considered longitudinal displacements might be different from the exactly 
solvable transverse case. 

The behaviour oi g{x) is put in perspective by the results for the distribution of near- 
est distances between particles and pinning centres shown in Fig.5. Either side of the sharp 
central peak there is a dip in the distribution showing the existence of a 'forbidden region' 
were the particles are less likely to sit. These regions disappear in the high temperature 
phase (i.e. when It > N) were the distribution is slightly peaked at the centre and decreases 
monotonically as r increases. We can account for the forbidden regions by considering the 
following simple model. We describe a single particle in the chain and the pinning centre 
nearest to it by the following model Hamiltonian 

H='^{x-xoy + Ap{{-^y-l) \x\<Rp 



.RpJ "-J l-^l ^ -'^P 



''■<=// r^ ^\2 



H=^{x-xoy \x\>Rp 

Thus we consider the particle to be interacting with a potential well of depth Ap, centred 
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at X = and of range Rp, and that this particle is coupled harmonically to the rest of the 
chain with an effective coupling constant fcg// and an elastic energy zero at position xq. 
The system is restricted in size to —a < {x, xq) < a, where the parameter a is given by the 
density of pinning centres in the simulations 2a = —. Since for each particle the parameter 
Xo will differ, when we calculate the distribution from the simulations, this is effectively the 
same as averaging the above Hamiltonian (or the probability function e~^^) over xq. Thus, 
we want to calculate 

I ra p-l3H{x,xo) 

where Z{xo,f3) is the usual canonical partition function with the exception that the x 
integration is over the finite interval —a < x < a. As usual (3 is the inverse temperature. 
The calculation reduces to a sum of error functions and the Xq integration has to be done 
numerically. The effective spring constant will be k^/f = ak, where l/a gives the number 



of springs which are affected by the displacement of a single particle |21, 22 



Fig. 5. shows the distribution P{r) obtained from a simulation with parameters 
N = 10^ Ap = 0.01, Rp = 0.1, Up = 1 and (3 = 10^. The inset to Fig. 5. shows the 
calculated distribution (P(x))^q with the same parameters as the simulation and with a 
renormalisation constant a = 0.175. The forbidden regions can be seen clearly in both 
cases. It is captivating to notice that the length scale 1/a = 5.7 is identical to the length 
scale Xc at which g{x) crosses between the two different algebraic regimes for this set of 
parameters. See Fig. 3 where the arrow at the x-axis locates the value log;^Q(l/a). 

It is of interest to measure the relaxation times of the system in order to decide upon 



the possibility of ergodicity, or replica, breaking in the disordered chain |T^. We measured 
the time dependence of the elastic energy Eei{t) (given by the first term in the Hamiltonian 
in Eq. ^. The relaxation of Eei{t) for different chain lengths is exhibited in Fig. 6 for 
Ap = 0.01, Rp = 0.1 at a temperature of T = 10~^. An average over different realisations of 
the random potential has been performed. The temporal behaviour of Eei{t) is independent 
of the length of the chain. Hence, we conclude that the disorder is unable to break the 
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ergodicity of the system. 



4 Conclusion 

We have studied a harmonic chain at low temperature subject to a random potential Fluc- 
tuations, both thermal as well as disorder, are very significant. This puts unpleasant limits 
on the system sizes which one can simulate in a reliable way. In spite of these limitations 
we conclude the following concerning the spatial behaviour (at low temperatures) of the 
displacement function g{x). A crossover length Xc exists. For x < Xc we find g{x) ~ x'^^'^ 
with rji ~ 3/4, in disagreement with the expected Larkin exponent. For x > Xc we find 
g{x) ~ x'^^^ where r]2 — 3/5, in agreement with replica theory. We are not able to iden- 
tify the functional dependence of Xc on the strength and range of the randomly positioned 
pinning wells. We do not find any difference in the behaviour of the open and the periodic 
chain. We are also unable to consider chains of sufficient length so to study the appearance 
of a third scaling region for g{x) > Cq. 

The average square fiuctuations W{L) ~ L^-^ was found to depend on the length of 
the chain with an exponent x consistent with the exact result for a transverse chain x = 2/3. 
The exponent might change on length scales beyond our reach. 

We studied the distribution P{r) of nearest distance between particles and pinning 
wells. As the temperature is lowered (or Ap is increased) particles become prohibited from 
a region around Rp away from the centre of the pinning wells. 

The crossover distance Xc and the distribution P{r) can be connected. In an effective 
single particle model of P{r) the scale Xc reappears as the ratio between the effective spring 
constant and the bare spring constant. 

The system was found to be self- averaging. A study of relaxation times in the system 
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showed that (as expected in one dimension) no ergodicity breaking occur. 
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Captions 

Figure 1. 

The displacement correlation function. Solid line corresponds to an average over 600 re- 
alisations of the random potential. The dotted line is generated by an average over 2300 
different initial starting configurations all annealed down into the same realisation of the 
random potential. Parameters are A^ = 200, Ap = 0.01, Rp = 0.1, np = 1, Tj = 1, and 

/ 



Tf = 10-3. 



Figure 2. 

A host of the concrete realisations oi g{x) entering into the averaged g{x) shown as the solid 
line in Fig. 1. 

Figure 3. 

Double logarithmic plot of g{x) in four different cases. Solid line represents Ap = 0.01 and 
Rp = 0.1. Dashed line corresponds to Ap = 0.01 and Rp = 0.2. The dot-dashed line belongs 
to Ap = 0.05 and Rp = 0.1. These three curves are for periodic chains. The dotted curve 
represents the measurement on an open chain with Ap = 0.01 and -Rp = 0.1. In all cases 
A^ = 10^, Tj = 1, and Tf = 10^^. The arrow indicates the location of log^gl^c) which 
coincides with the location of log^glV'^)- See the text for an explanation. 

Figure 4. 

Double logarithmic plot of the average fluctuations in the displacement W{L) as function 
of L (crosses) together with the amplitude g{L/2) of the displacement correlation function 
(circles). 

Figure 5. 

The distribution of nearest distance between particles in the chain and pinning centres. 
Parameters are Ap = 0.01 for the solid line and Ap = 0.05 for the dashed curve. In both 
cases Rp = 0.1, Up = 1, N = 10^, Tj = 1, and Tf = lO^'^. The insert corresponds to the solid 
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curve. It is calculated from the single particle model with a scaling paramenter a = 0.175. 

Figure 6. 

The time dependence of the elastic energy (measured per particle). The simulation is run 
at fixed temperature T = 10^^. Curves for A^ = 200, 300, 400, and 800 are shown. In all 
cases Ap = 0.01, Rp = 0.1, Up = 1. 
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